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Abstract. We present a linear analysis of the stability of accretion disks in which angular momentum is removed by the 
magnetic torque exerted by a centrifugally driven wind. The effects of the dependence of the wind torque on field strength 
and inclination, the sub-Keplerian rotation due to magnetic forces, and the compression of the disk by the field are included. 
A WKB dispersion relation is derived for the stability problem. We find that the disk is always unstable if the wind torque is 
strong. At lower wind torques instability also occurs provided the rotation is close to Keplerian. The growth time scale of the 
instability can be as short as the orbital time scale. The instability is mainly the result of the sensitivity of the mass flux to 
changes in the inclination of the field at the disk surface. Magnetic diffusion in the disk stabilizes if the wind torque is small. 
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1. Introduction 

Magnetically driven winds from accretion disks have been 
considered as promising explanations for highly collimated 
jets observed in many classes of astrophysical objects, rang- 
ing from active galactic nuclei (AGNs) to young stellar ob- 
jects (YSOs)(see the reviews by Blandford 1993, 2000; Pringle 
1993; Spruit 1996). Hydromagnetic winds from accretion disks 
provide a mechanism to produce collimated jets as well as an 
efficient angular momentum loss mechanism from accretion 
disks. The jets are powered by the gravitational energy released 
by accretion of matter through a magnetic field with an open 
configuration. If the inclination angle of the magnetic field lines 
is sufficient large with respect to the axis of the accretion disk, a 
magnetically controlled, centrifugally accelerated wind is nat- 
urally driven from the disk surface. Open magnetic field lines 
threading the accretion disks are a crucial ingredient in this 
model (Blandford & Payne 1982). A global, ordered magnetic 
field is usually assumed in investigations, though the origin of 
the ordered magnetic field that is assumed to thread the disk is 
still not clear. One possibility is that the field is advected in- 
wards by the accreting matter in the disk (Blandford & Payne 
1982). The other one is that the field is generated locally by a 
dynamo (e.g. Tout & Pringle 1996). A steady open magnetic 
field may be maintained if the inward dragging of field lines 
by the disk is balanced everywhere by the outward motion of 
field lines due to magnetic diffusivity (Lubow, Papaloizou, & 
Pringle 1994a). Lubow, Papaloizou, & Pringle (1994b, here- 
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after LPP) explored the stability properties of the disk-wind 
system on the assumption that the angular momentum of accre- 
tion disk is removed by the magnetic torque alone. They argue 
that such a disk wind is unstable. Konigl & Wardle (1996) on 
the other hand claim stability of their own disk-wind solutions 
(Wardle & Konigl 1993). Krasnopolski et al. (1999) also claim 
stability on the basis of numerical simulations. Instability of the 
type proposed by LPP has been incorporated into a numerical 
model by Agapitou and Papaloizou (in preparation). 

2. Model 

Magnetically driven winds can remove angular momentum 
from the disk. The wind acceleration is sensitive to the inclina- 
tion of magnetic field at the disk surface (Blandford & Payne 
1982). If the magnetic field is advected inwards by the accreted 
matter, the field inclination at the disk surface is mainly de- 
termined by radial velocity of the matter in the disk. A small 
perturbation increasing the radial velocity causes the poloidal 
field to be bent closer to the disk surface, so that it gives rise to a 
higher mass loss. By the increased angular momentum loss, the 
radial velocity continues to increase. Thus an instability arises 
(LPP). 

A number of factors have to be included in the analysis. If 
the magnetic field strength increases inward in the disk, it ex- 
erts an outward radial force providing partial support against 
gravity, so that the rotation of accretion flow is sub-Keplerian. 
The height of the disk is compressed by both the vertical com- 
ponent of gravitational force and the vertical pressure of the 
radial magnetic field component. 
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In principle, one can numerically calculate the disk struc- 
ture and the magnetic field configuration in the disk simultane- 
ously (Ogilvie 1997, Ogilvie & Livio 1998). For simplicity, in 
this work we use an approximate field configuration to estimate 
the scale-height of an isothermal disk compressed by the field. 

We assume that the accreting matter is corotating with mag- 
netic field lines in the disk. This is an appropriate limit, if the 
field is sufficiently strong (e.g. Ogilvie & Livio 2001). The mat- 
ter at the disk surface is centrifugally accelerated along the 
magnetic field line anchored at the disk surface. The matter 
corotates with the magnetic field line roughly till it approaches 
Alfven point. Beyond the Alfven radius, the magnetic field ex- 
erts no further torque on the escaping wind matter. The mag- 
netic torque can be estimated if the mass loss rate in the wind 
and the Alfven radius are known. As usual, the mass loss rate 
in the wind is determined by the density at the slow sonic point 
in the wind. The slow sonic speed is close to the sound speed, if 
the Alfven radius ta 3> r s . To avoid solving the equation of ra- 
diative transfer in vertical direction of the disk, we approximate 
the temperature of the disk to be a constant along the vertical 
direction. 

For most of the stability analysis, a local approximation 
can be made involving the values of physical quantities near 
the foot point of a field line. The analysis, however, involves 
the magnetic torque in an essential way, hence the value of the 
Alfven radius r A for the field line must be known. Since the 
Alfven radius depends on the global field configuration, not 
only on the field strength at the base of the field line, it can be 
regarded as an external parameter for the stability analysis. To 
minimize the number of additional parameters in the analysis, 
we assume the field at the Alfven surface to be roughly self- 
similar, varying as a power of distance from the axis. The value 
of rA can then be determined and we can evaluate the torque 
exerted by the wind on the disk. It depends on the mass flux 
along the field line, which depends on the density at the sonic 
point. This in turn depends sensitively on both the temperature 
of the disk and on the inclination and strength of the field near 
the disk. 

Changes in torque exterted cause the mass distribution in 
the disk to change. This in turn causes changes in the distribu- 
tion of the field lines anchored in the disk. Since the mass and 
angular momentum loss rate on a field line depends sensitively 
on the field configuration (in particular the inclination), there 
is a strong intrinsic feedback, and it is the purpose of our in- 
vestigation to see under which conditions this feedback leads 
to instability of the type proposed by LPP. 

The stability problem is in principle global in nature. We lo- 
calize it by a WKB approximation (short wavelength limit) and 
a dispersion relation is then derived for the stability problem. 
This is analogous to the so-called tight-winding approximation 
used in studies of self-gravitating disks. 

A further simplification we make is that the magnetic 
torque in a steady wind is also be used in the perturbations. 
This assumption is valid if the time scale of the instability is 
much longer than the disk-wind coupling time scale. We dis- 
cuss this further in the last section. 

Since our focus is on the instability caused by the wind 
torque, we leave out angular momentum transported by the 



usual disk viscosity, so that viscous modes are also absent in 
stability analyses. The radial pressure gradient in the momen- 
tum equation is neglected for similar reasons, to avoid the pres- 
ence of an acoustic mode. The circular motion of the disk in 
equilibrium is given by the balance between the gravity and ra- 
dial magnetic force. To evaluate the perturbations in the mag- 
netic field above the disk we assume this field to be potential. 
This a good approximation as long as the Alfven surface is far 
from the disk, since the short-wave approximation ensures that 
the magnetic perturbations decay quickly with height above the 
disk. The perturbed magnetic field configuration above the disk 
is then determined by the normal component of the perturbed 
field at the disk surface. This, in turn, follows from the radial 
displacements in the disk, through the induction equation. With 
respect to the magnetic field and the forces it exerts, the model 
is thus the same as that used earlier by Tagger et al. (1990), 
Spruit and Taam (1990), Lubow and Spruit (1995), Spruit et al. 
(1995), and Stehle and Spruit (2001). 

To evaluate the mass loss rate in the magnetic wind, a 
model for the distribution of mass with distance above the mid- 
plane is needed. We assume the Gaussian distribution valid for 
an isothermal disk, but also include the effect on the disk thick- 
ness of vertical confinement by the magnetic pressures. 

Strictly speaking, the presence of a magnetic torque will 
lead to an azimuthal component of the field. However, the az- 
imuthal field component at the disk surface is always small 
compared with the radial component for any wind models as 
long as the Alfven radius is far from the radius of the field foot- 
point at the disk mid-plane (see further discussion in section 
3.2). If the magnetic field strength inside the disk is suffiently 
weak, differential rotation acting on the r- and z-components 
will produce a time-dependent azimuthal field. Such a nonsta- 
tionary state will lead to magnetic turbulence such as found 
in numerical simulations of an initially weak field (Hawley et 
al. 1995). The field we consider here, on the other hand, is 
a systematic poloidal field (i.e. crossing the disk plane), and 
sufficiently strong such that Balbus-Hawley instability is sup- 
pressed. At this strength, the field is also sufficiently strong 
to enforce 'isorotation' : constancy of the rotation rate along 
field lines. Though the origin of such fields is not as clear as 
the turbulent fields that evolve from initially weak seed fields, 
they are by far the most logical configurations for producing 
systematic outflows from accretion disks. Spruit et al. (1995), 
Stehle (1997) and Stehle and Spruit (2001) have shown that 
even surprisingly strong poloidal fields, approaching (a frac- 
tion of) equipartition with orbital kinetic energy, can still be 
stably anchored in a disk. These analyses did not include insta- 
bilities associated with magnetically driven outflows however, 
which are the subject of the present study. 

3. Magnetic torque 

From standard magnetic wind theory, we know that the angular 
momentum flux along a field line is given by T — rh£l(r{)r\, 
where r is the cylindrical radial coordinate, ta the Alfven ra- 
dius, m the mass loss rate, and Q the rotation rate of the foot- 
point at radius n on the disk. This is the total angular mo- 
mentum flux, including both the magnetic torque and the flux 
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carried by the mass itself. For the effect on the disk, only the 
magnetic torque T m is consequential; evaluating it at the disk 
surface, and counting both surfaces, we have 



2rMVn)(ri-rf), 



(1) 



In the limit rA > T m ~ 2rnQ.Q{r{)r\. 

The mass loss rate in the wind is governed by the position 
of the sonic point: 



B z 



(2) 



where z s is the height of the sonic point above the midplane of 
the disk. 

As discussed in section 2, the mass loss rate of wind from 
an isothermal disk is estimated as 



m • 



) poc s exp[-(* es -* ei )/c s 2 ], 



(3) 



where ^cs and 'J'ei are the effective potential at the sonic point 
and the footpoint of the magnetic field line respectively. The 
effective potential \& e is given by 



*e(r,«) 



GM 



1 



( r 2 +z 2)l/2 2 



Q 2 (n)r 2 



(4) 



where r ; is radius of the magnetic field line footpoint at mid- 
plane of the disk. 

The position of the sonic point can be estimated, if we know 
the shape of the magnetic field lines between the footpoint 
and the sonic point. In principle, the field line shape is com- 
putable by solving the radial and vertical momentum equations 
together with suitable boundary conditions. Analytic solutions 
are available only for two extreme cases, the weak and strong 
field extremes. We find that the field line shapes are slightly 
different in these two extreme cases. In the strong field case, 
the disk is compressed mainly by the vertical component of 
magnetic force, and the shape is similar to the Kippenhahn- 
Schluter model (KS) (Kippenhahn & Schluter 1957), valid for 
an isothermal sheet suspended against gravity by a magnetic 
field. The exact shape of the field lines depends on the ther- 
mal structure of the disk. For our purposes an isothermal layer 
would be a sufficient model, and the KS model applicable. Its 
analytic form is still a bit too cumbersome however, so we have 
approximated it by the following expression: 



r — n = 



H 



H 



(1 - V f + rg*H-*YI* - -^(1 - ?7i 2 ) 1/2 ,(5) 



«0»?i 



where H is the scale height of the disk, j]\ = tanh(l). The in- 
clination of the field line n — B z /Bf is a function of z. By the 
index o we denote the value kq measured at the nominal disk 
surface z = H. High above the disk, the inclination then be- 
comes kq tanh(l). We find that this expression reproduces the 
basic features of the KS model well for both weak and strong 
field cases. 

We approximate the position of the sonic point as the max- 
imum of the effective potential along the magnetic field line. 
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Fig. 1. The magnetic torque as functions of the dimensionless 
angular velocity of the disk and the magnetic inclination k , 
for a disk thickness Hj r = 0.01. The angular velocity is also a 
measure of the field strength (1 - (£1 /^k) 2 ~ B 2 ). 



The height above the disk of the sonic point at radius r of the 
disk is then from using the potential (4) and the magnetic field 
line shape (5): 



H 



{ 



3n 2 H - k 2 t] 2 H 



2U/2 



1/2 



l + »7f 



,(6) 



which is valid for z s /ro <C 1. 

To complete the estimate, we need the value of the scale 
height of the disk H , which can be calculated from the vertical 
component of momentum equation 



; dp{z) 
dz 



B r {z) dB r {z) 



4.7T 



dz 



(7) 



The density distribution p(z) in z-direction can be obtained by 
integrating Eq. (7) with the given field configuration and the 
assumed disk temperature. 

As a result of the magnetic forces the distribution is not ex- 
actly the Gaussian distribution of an isothermal disk. The disk 
thickness is reduced by the magnetic forces. In this case, we 
define the disk scale height H as p{H) = p(0) exp(— 1/2), the 
disk scale height H can then be evaluated numerically. From 
these numerical values, we have constructed a fitting formula 
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where B z (r n ) is the z-component of the magnetic field at the 
disk surface. 

At the Alfven point, the Alfven velocity is 



(47rp A )V2< 



(11) 



As in all magnetic acceleration models, its value is of the order 
V A ~ r A^o, where B^ and p\ are the poloidal magnetic field 
strength and density of the outflow at Alfven point. 

Mass and magnetic flux conservation along a magnetic field 
line requires 



777 



(12) 



Combining Eqs. (10)-(12), the Alfven radius is found: 



47T777£!o(7'o) 



(13) 



Fig. 2. Same as Fig. 1, but for H/r = 0.001. 

used in the rest of the calculations to represent the parameter 
dependence of the scale height H: 



H=- 



4(l-C- 1 /2)2 K 2^4 



1/2 



4(1 -e-V2)«of22 * 



(8) 



We find that the scale height H given by expression (8) is ac- 
curate to about 1 % compared with the numerical results in the 
range of the parameters adopted in this paper. The actual den- 
sity distribution deviates from a Gaussian distribution, so the 
density at the sonic point still is only approximate, but the fit- 
ting formula properly takes into account the compression of 
the disk in the strong field case. Note that Eq. (8) reduces to 
H = c s /£Ik in the Keplerian case. Now, using Eqs. (4)-(6), we 
obtain the effective potential difference A^ c between the sonic 
point and the magnetic footpoint in the mid-plane of the disk 



A* = „ 



1 



O 2 7 2 



+ 



(g| - Sl 2 Q )r H 



(i-tf + tf#r-»)V a -(i-»J?) 1/a 



H 2 

2 4 



Now we can evaluate the magnetic torque T m by using Eqs 
(3), (6) and (13). 

We define dimensionless quantities by 



(l-rg + rg^H-^-il-ri) 1 



/2 



Using Eqs. (3)-(9), we can evaluate the mass loss rate in wind. 
The magnetic field far from the disk surface is assumed to be 
roughly self-similar as discussed in section 2: 



-Bp ~ B z (r )f 



B z (r ) 



,a > 1, 



(10) 



* = h = —. 



m S ^r 2 ' 
The dimensionless magnetic torque T m is then given by 

f m = o|^(i - o 2 )^tt k^ 1 k°*x (i + K 2 s y&y 



(14) 



2(1 - e-V 2 )KoH 



2(a + l) 



a- I A* e \ 

<xp 1 ' ' (15) 



The parameters of the model are now H, Cl , k and the self- 
similar index a of magnetic field shape. Here k s is the inclina- 
tion of the field line at the sonic point. 

3. 1. Dependence of the magnetic torque on field 
strength and inclination 

With increasing field strength, the rotation of the disk becomes 
more sub-Keplerian. Though the increasing field strength it- 
self increases the torque through the increasing Alfven radius 
fse^nd factor in Eq. (15)], the lower rotation of the footpoints 
also increases the potential barrier that the wind has to over- 
come (last factor). In many cases of interest, an increasing 
field strength will reduce the magnetic torque (e.g. Shu 1991; 
Ogilvie & Livio 1998, 2001). In the following, we use the ro- 
tation rate relative to Keplerian, f2 , as measure of the field 
strength. Figures 1 and 2 show expression (15) for the cases 
H = 0.01 andff = 10~ 3 . 
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3.2. Azimuthal magnetic field component at the disk 
surface 

The magnetic torque T m exerted on the unit surface area of the 
disk is given by 



T — 



B?Blr 



2tt 



where B^ is the azimuthal field component at the disk sur- 
face. Let g m be the radial acceleration due to magnetic ten- 
sion, averaged over the surface mass density of the disk, g m = 
BfB z /(2irT,). The radial balance of forces is then rfi 2 , = 
g — g m . Hence we have the relation 



3 
B* 



(17) 



For the parameter values where our analysis will be valid, the 
azimuthal field strength at the disk surface is small, B^ <C Bf 
(see Figs. 1 and 2). 

4. Disk equations 

As discussed in section 2, we neglect the usual disk viscosity 
in the problem. Torques due to the magnetic wind are therefore 
the only source of angular momentum loss of the disk material, 
and in their absence the mass flux through the disk vanishes. 
The mass conservation and angular momentum equations for 
the disk are then 



— +v —\ £ + £— 

dt r dr I dr 



and 



d_ 

dt 



dr 



V<f, 



T 

- 1 - m 



(19) 



where T m is the torque contributed by the magnetically driven 
wind per unit area of disk surface. The radial component of the 
momentum equation is 



d_ d_ 

dt r dr 



v r + g 



r 



0. 



(20) 



where the gravitational and radial magnetic forces are given 
by g = GM/r 2 and g m = BfB z /(2nT,), and we neglect the 
gas pressure gradient as discussed in section 2. The magnetic 
induction equation is 



d d \ v r B z dv r n d ( dB r , 
- + v r —\ B z +^^+B z -^ + ^— r—^ = 0.(21) 
dt dr I r dr r dr \ dz 



where r\ is magnetic diffusivity, H is the disk scale-height. 

As discussed in section 2, a potential field above the disk is 
assumed, so that the stream function of the field satisfies 



<9 2 $ 
dr 2 



1 d<P 

r dr 



9 2 $ 
W 2 



= 0. 



(22) 



4. 1. Linearized equation in a local approximation 

As discussed in section 2, the problem can be localized by a 
WKB approximation in a short wavelength limit. The perturbed 
stream function is still given by Eq. (22): 



<9 2 <5$ 1 d5<S> d 2 5<$> 

+ -^r- = 0. 



(16) dr 2 r dr 



dz 2 



(23) 



The strength B z (r) at the disk surface acts as boundary condi- 
tion for this potential problem. Its solution is in general a global 
one, the solution at any point (r, z) depends on the entire field 
strength distribution B z (r). In a short wavelength limit, how- 
ever, the problem becomes local again. Our short wavelength 
limit is the magnetic equivalent of the so-called tight-winding 
limit in self-gravitating disks. 

Equation (23) has separable solutions 5$>(r, z) — 
F(r) exp(— k\z\), where F(r) satisfies 



F" - -F' + k 2 F = 0. 

r 



(24) 



Let y = k(r — r ). In the short-wave limit kr ^> 1, we have 

(25) 



d 2 F dF p _ {) 
dy 2 dy 



where e — 1/kr (e <s£ 1). 

Denote Eulerian perturbations of a quantity q by Sq, 
Lagrangian perturbations by an index i, and the equilibrium 
state by an index - 

The perturbed fraction of the stream function is: 



(18) 5®(r,z) = Acxp ±iy - k\z\j . 



The perturbed components of the magnetic field are 

6B r (r,z) = -^ = -6*(r,z), 
r dz r 

and 

SBz(r,z) = - - ( - +ik \ 5*{r,z). 



(26) 



(27) 



(28) 



Here we have taken the sign in Eq. (26) to be "+". The mag- 
netic field components are related by 



5B Z = 5B r iSB r- 



(29) 



Let a dot be the Lagrangian time derivative, and £ be the 
Lagrangian displacement. Neglecting the gradient of r in the 
background and assuming H to be a constant along r, the lin- 
earized equations of the disk are: 



El + H drZ T - 0, 

& + (2fi + S)£ + ^-SgE 1 = 
where the shear rate S = r{dQ /dr). 

£ r - 2£! O £0 - 5ml = 0, 



(30) 
(31) 

(32) 
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where rSl 2 , = g - g m , 



B Z 1 + B z ()drir + -JjB z \ 



0. 



(33) 



As discussed in section 2, we assume that the magnetic 
torque as evaluated for steady conditions are also valid for the 
perturbations. The magnetic torque T m is a function of the an- 
gular velocity of the flow £1 and the magnetic field inclination 
Ko, «o = B z /B^, at the disk surface. The perturbed fraction of 
the magnetic torque ST m is given by 



STm = dT in ( K ,n) SK+d T in ( : ,n) §n 



dn dfl 
dT m (n ,n) ( 5B Z t ,B z0 6B z \ t oT m (/t,Q) 



B r n B 2 



rO 



sn 



on 

. 2 5B Z dTj^t) 

d. ^ + ik ° ) bT + ^^ sq - 

where Eq. (29) is used. The perturbed fraction of the radial 
magnetic force is given by 



5. Dispersion relation 

We write the perturbed quantities as 

£ r ~ exp[i(u;t + kr)]. 

The dispersion relation follows from Eq. (37): 



4 . ( I <9T m rjk\ 3 



Koffmofc + 2f2 (2£!o + S) 



Z r 2 H dfl 



2kn dT m , t 2 : .gmoKokdTn- 



-(kq + in ) + i- 



S r dn v uy £ r 2 dn 

2kfl T m0 T]k 2 g ma 2rjkn 



+ 



S r 

_ . yk 2 g m0 dT m 2r]k 2 noT mQ 



= 0. 



(38) 



'5 B r SB Z 5E 

0g m = SmO I + d 7T~ 



(34) 



Neglecting the radial gradients of the equilibrium quantities 
(a local approximation) the Eulerian and Lagrangian variations 
are equivalent, so that the perturbed torque and magnetic accel- 
eration are 

T m i = « («o + ««o)^— + ^ (35) 

and 

5ml = 5m0 (1 - IK )-^- + <9 r £ r , (36) 



where Eq. (30) has been substituted. Combining Eqs. (30)-(33), 
we have 



+ 



£ r 2 dCt H , 



2f7 (2fi + 5) 



We use some dimensionless quantities in addition to those de- 
fined in Eq. (14): 

A-— ill—— S - r ° 

H ' Six ' dr 



ffmO 



9m0 



and the quantities p, q: 



V 



(39) 



dT m _ p ~ 

~ ft" m0 ' 



t: — — — J m0, 

OK Kq 



which measure the dependence of the wind torque on the rota- 
tion rate and field inclination, respectively. 
Equation (38) becomes 



. ( p - 2tt?7\ 



^(l-n§) 

XH 



+2Sl (2O + S) + 



2irpfjT ma 
tin 



u 2 



4irQ T m0 
XH 



20o dT m . . 2 . 2iloT m o .SmoKo STjx 

(Kq + «K ) - h «- 



+ 



S r 



S r 2 dn 



H 



d r i r + ^^(2n + s)i r 



2irpK T m0 ~ 2 4irqn 

-1 -- - — (1 - Sin) 

ASl ff XH 



T m0 (l + *Ko) 



A 2 i/ V 07 A 



' vkgr^dTyn _ 2r)kn T m0 
S r 2 H 90 S riJ 



9 r ^ r = 0. 



(37) / '" J ' // ' /|M -:l L' 
A 2 il i? 



, 1 A2^ 1 „.87r 2 77il T m o 



A 2 i? 



= 0. 



(40) 
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5. 1. Restrictions on parameter values 

The dispersion relation contains several small parameters, since 
a number of conditions have to be satisfied for consistence with 
the assumptions made. For a local approximation to be valid, 
we must have 



A/r < 1. 



(41) 



On the other hand, since we have approximated the disk as thin, 
the wavelength has to be large compared with the scale height: 



5.2.1. fj ^0,T m0 = 

We now analyze the marginal stability, i.e., a); = 0. In the ab- 
sence of the magnetic torque (T m0 = 0), Eqs. (43) and (44) 
reduce to 



3tf 



-^^i + 2Sl (2ilo 
A 



+^(i 

XH 



A 2 iJ 



Q 2 ) = 



(45) 



H < A. 



(42) 



An unknown in the problem is the magnetic diffusivity, 
hence we must explore the dependence of the results on 77. In 
our assumptions, we have ignored viscosity in the disk. If the 
magnetic diffusivity is due to the same small scale process that 
causes the viscosity, it is then most consistent to set 77 = as 
well. On the other hand, it is instructive to explore the conse- 
quences of allowing some diffusion, since it is likely to damp 
instabilities. With our scalings (39), the relevant value is then 
fj ~ a ~ 1, where a is the viscosity parameter. We report 
results for both cases, fj <C 1 and fj ~ 1. 

5.2. Approximate analytic results 

The nature of the unstable modes can be explored by looking 
at limiting cases. We can look at the case where the magnetic 
torque is weak, T m0 <C 1. Assuming 1/A and T m0 to be the 
same order and ignoring quadratic and higher terms, the last 
two terms in Eq. (40) containing T m0 / A can be omitted. The 
dispersion relation reduces to a third-order equation. Let u> = 
u) T + iu>i, it can then be separated into two equations: 



3u>f — 2u>[ 



, fi A 



+ 



2irpfjT m0 
fin 



XH 



£>r "I ~ g - " ^mO 

XH 



(43) 



and 



27T77 _ 9 



3w 2 



2fi n (2fin + S) 



-^(1-^) 
XH 



Qi + 2 -^tf-^(2n + S) = 0.(46) 
A A 



Let Wi = in Eqs. (45) and (46), we have 

w r = ±[2O (2f) + £)] 1/2 , 

and 



V : 



^oA . 



(47) 



(48) 



Since 77 > 0, cD r /A has the same sign as kq. Since the field 
strength will normally (but not necessarily) decrease outward 
in the disk, the inclination of the field lines will be away from 
the axis, i.e. k > 0. In the absence of diffusion, a transition 
from stability to instability, if it exists at all, then is possible 
only for an inward traveling wave mode, i.e. u> r /A > 0. 

Further progress can be made if we assume that the dif- 
fusivity is small, and by looking at conditions near marginal 
stability. Treating fj it as a small quantity of the same order as 
the other small quantities, i.e. 

fj ~ Cbi ~ 1/A, 

the growth rate is then given by 



XH(3Cu? - wg) ' 



(49) 



where 



and 



pT^ + 2^ 



fin 



3^ r 2 - 2fi (2fi + S) - ^0 (1 - fi 2 ) 2 -^I^ 
XH fi 



ff^(l -fi 2 )f m0 -l*(2fi + S) = 0. 
XH\Iq X 



(44) 



= 2fi (2fi + S) + 27tko(1 - til)/XH, 
A 1 = 2wfjH^ - 4Trf]Hfl a (2n + S). 

The real parts of the roots of the equation are: 
Ci 



A 2 £ 2 tf' 



where 

Ci = 4tt 2 ?7(1 - 
and 

<D r 2,3 - ±^0 - 



oh 



2X 2 HCj 2 



(50) 
(51) 

(52) 

(53) 
(54) 
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The imaginary part of u> for these roots is 



and 



Wil ~ 



and 



Wi2,3 = 



47rrySl (20o + S') 



Aw 2 



2X 2 HQ 2 ' 



(55) ^r2,3 - ±^0 - 



C 2 



(56) 



The modes 2, 3 have the same growth rate in the limit taken 
here, where only the first order in fj is kept. From Eqs. (55) and 
(56), we find that all three modes are stable in low-?7 limit. It 
implies that the neutral wave and the outward traveling wave 
are always stable. 

The stability condition for the inward traveling wave mode 

is 



ry < ^[2ft (2f)o + S)] 



1/2 



(57) 



Since the inclination n is of order unity, and fj,l/X are small 
numbers, the stability condition is satisfied. Not surprisingly, 
the disk is stable in the absence of a wind torque. 

5.2.2. fj = 0, T m0 / 

We now include the magnetic torque. In this section we treat it 
as a small quantity of same order as the other small quantities in 
the problem. In the absence of the magnetic diffusion {fj = 0), 
Eqs. (43) and (44) reduce to 



3tf - + 2 f2 (2f2 + S) 



XH 

and 

- 3 pT m - 2 

UJ: = LU: - 

O 



+ ~~ -Tmo ~-~-T m0 = (58) 

XH XH 



3J> r 2 - 2f2 (2f2 + s)- ^^(l - o 2 ,) 

XH 



+ — =; — w r — J m0 - --- (1 - S2 )i m0 - U.(59) 

sZq A-rz AHUq 

We analyze the situation close to marginal stability in the 
case of low magnetic torque (T m0 — > 0). The growth rate is 
given by 



A 3 



(60) 



Wl XM (3Cu2 -wg)' 
where 

A 2 = ~XH P f mQ u 2 r - 47rg/t ^o^mO - 27r P K (l - ng)T m0 .(61) 
The real parts of the roots of the equation are: 

C 2 



where 

C 2 = 47rqfi f , m0 - 47rOoT , m0 



(62) 



(63) 



2XHCjI 



The imaginary part of <I> for these roots is 

£>ii = [4^g Ko ^^ + 27rp Ko (l - ng)] J^ 



and 



^i2,3 = [pAff (2ft + S) - 2^ Ko ^o]-^. 



(64) 



(65) 



(66) 



The modes 2, 3 have the same growth rate in the limit taken 
here, where only the first order in 1/ A and T m0 is kept. 

5.3. Interpretation of the analytic results 

The growth rates of the three modes (65, 66) are all propor- 
tional to the magnetic torque. Hence they all qualify as instabil- 
ities caused by the coupling between the disk and the outflow, 
if their Wj are negative. The first mode, u>\ is neutral, in the ab- 
sence of diffusion and magnetic torques. If only magnetic dif- 
fusion is included (section 5.2.1), it is a damped mode. Thus, 
low frequency perturbations are stabilized by magnetic diffu- 
sion, in the absence of a magnetic torque. When a magnetic 
torque is present (section 5.2.2) the sign of the growth rate is 
determined [see eq (65)] by the quantity 2q£l 2 + p(l — fl 2 ). 
For low field strength, the second term is small, and instability 
is determined by the sign of q, the rate of change of the mag- 
netic torque with inclination of the field line. As Fig. 1 shows, 
the magnetic torque increases with decreasing inclination (field 
lines bent further away from the vertical). This reflects the fact 
that the centrifugal effect is stronger on field lines that are bent 
further away from the vertical, increasing the mass flux and 
magnetic torque. Thus q is negative in the cases of interest, and 
the growth rate — u)\ positive. The term involving p, the depen- 
dence of the wind torque on rotation rate, is stabilizing since 
p > (cf Fig.l). In dimensional terms, the approximate condi- 
tion for instability in the absence of magnetic diffusion and in 
the limit of a weak magnetic wind torque is thus: 



-dT u 
8k 



> 



9r 



(67) 



where g m — (f^ — fl 2 )ro = B z j (27r£) is the acceleration 
due to magnetic support of the disk against gravity, (the minus 
signs are included because the derivatives are negative in the 
cases of interest). 

This mode evidently represents the sought instability, as- 
sociated with the dependence of magnetic torque on field line 
inclination. The result shows, however, that the instability is not 
governed only by the dependence of the torque on inclination. 
If the magnetic field is strong enough to change the rotation 
rate by providing support against gravity, the second term in 
(65) is important, and its effect stabilizing. This represents the 
fact that a lower rotation rate increases the height of the po- 
tential barrier that the wind has to overcome (cf discussion in 
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section 3), which decreases the magnetic torque (Ogilvie and 
Livio, 2001). 

We briefly discuss the two oscillatory modes. Since p > 
and q < for the cases of interest (Fig.l), eq. eq. (66) shows 
that these modes are stable. In the limit of i) = and T m0 = 0, 
the two modes have the same frequency (cf. Eqs. (54) & (64)) 



Wr2,3 - ± 



2tt (2ft + S) + (1 



tt 2 



-,1/2 



(68) 



[2f2 (2^o + S)] 1 / 2 is the epicyclic frequency associated with 
the stable angular momentum gradient. The second term repre- 
sents the restoring force due to the change of inclination of the 
field lines by the perturbations (note that 1 - f) 2 , oc B r B z ). On 
its own, this magnetic restoring force would yield a wave with 
frequency uj oc A -1 / 2 , like a surface wave. This is the result of 
the long-range nature of the magnetic perturbations of the po- 
tential field outside the disk. These magnetic waves have been 
analyzed before (Spruit & Taam 1990; Tagger et al. 1990). 
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6. Numerical Results 

We obtain the growth rate of the instabilities by solving Eq. 
(40) numerically. For the shear rate S we take the Keplerian 
value S = -3/2. 

The problem is described by the dimensionless disk scale- 
height H, the angular velocity of accretion flow £7 which is a 
measure of the magnetic field strength, the inclination of mag- 
netic field k = B z /B r at the disk surface, and the dimension- 
less magnetic diffusivity rj. The inclination turns out to influ- 
ence the stability mainly through its effect on the wind torque, 
which depends sensitively on it (Figs.l, 2). The second main 
parameter determining the stability is the strength of the mag- 
netic field, which we have expressed throughout this paper by 
its effect on the rotation rate (relative to the Keplerian rate) and 
the angular velocity of accretion flow, f) = ft /fix- With in- 
creasing field strength, the magnetic stresses start contributing 
to the support of the disk against gravity, decreasing the rota- 
tion rate. 

We take the self-similar index a of the magnetic field shape 
to be 4 in the following. In Figs. 1 and 2, the magnetic torque as 
functions of k and f2 is plotted for different values of the di- 
mensionless disk scale-height H. For small disk scale-height, 
i.e., if the temperature of gas in the accretion disk is low, the 
magnetic torque is significant only if the angular velocity is 
close to Keplerian velocity and the magnetic field inclination 
angle to the disk surface is low (low no, see Fig. 2). The nu- 
merically determined stability boundaries for the neutral wave 
mode are plotted in Figs. 3 and 4. It is seen that instability hap- 
pens in basically two region of parameter space. One is at high 
wind torque (right side of the diagrams); the instability then 
happens roughly independent of the strength of the magnetic 
field. At low wind torques, instability is possible only when the 
rotation rate is close enough to Keplerian. 

Examples of the growth rate of the neutral wave mode for 
different parameter values are shown in Figures 5-11. 



Fig. 3. The stability boundary of the neutral wave mode, for 
disk thickness H/r = 0.01 and wavelength X/H = 10, as 
a function of wind torque and deviation from Kepler rotation, 
the strength of the magnetic field through the radial equilibrium 
condition, with high field strength corresonding to low Q/Qk- 
The dimensionless diffusivity is i) = 0.75(solid line), 0.25(dot- 
ted) and O.Ol(dashed). 
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Fig. 4. Same as Fig. 3, but for H/r = 0.001. 
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Fig. 5. The growth rate of the instability as functions of mag- 
netic field inclination k and the dimensionless wavelength 
X/H. The dimensionless shearing rate S = —3/2, disk thick- 
ness H/r = 0.01, magnetic diffusivity ij = 0.75 and angular 
velocity fi /^K = 0.995 are adopted. 




Fig. 6. Same as Fig. 5, but for different value of angular veloc- 
ity: n /n K = 0.99. 




Fig. 7. Same as Fig. 5, but for n /^K = 0.9. 

7. Discussion of the results 

We have presented a linear stability analysis of disks with mag- 
netically driven winds, and confirm the instability by the mech- 
anism proposed by LPP. 



/• 5 r 




Fig. 8. Same as Fig. 5, but for fi /O K = 0.85. 
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Fig. 9. Same as Fig. 6, but for H/r = 0.001. 

In the presence of a strong magnetic field, the disk is com- 
pressed by the vertical component of gravitational force as well 
as vertical magnetic pressure, and the vertical structure of the 
disk is significantly altered, which affects the position of sonic 
point and the wind torque. In this strong field case, the angu- 
lar velocity of accretion flow significantly deviates from the 
Keplerian velocity due to the radial magnetic force, and the 
magnetic torque is then negligible for any magnetic field in- 
clination angle (see Figs. 1 and 2, see also Ogilvie and Livio 
1998). We have taken these effects into account in present in- 
vestigation. 

The system of equations analyzed has four modes. In the 
absence of a magnetic field, two are a neutral displacements. 
One of these two is not specifically related to the magnetic wind 
torque, and we have excluded it from quantitative analysis. The 
other has a frequency proportional to the magnetic torque (real 
and imaginary parts of the frequency are of the same order of 
magnitude), it can become unstable. 

The final two modes represent an inward and an outward 
traveling wave. The restoring force in the wave is a combina- 
tion of the coriolis force (epicyclic motion) and the magnetic 
forces. Both inward and outward traveling waves are stable in 
the range of validity of our assumptions. 

In Fig. 5, the growth rate of the instability is plotted for 
angular velocity very close to Keplerian velocity (Oo = 0.995), 
which may approximate the case considered in LPP. The disk 
becomes more unstable if the angular velocity of the flow is 
close to Keplerian (see Figs. 5-8). For high inclination angles 
of magnetic field line with respect to the surface of the disk 
(large Kq), the magnetic torque is very small. Instability is then 
suppressed by magnetic diffusion. 



Fig. 10. Same as Fig. 6, but without magnetic diffusion, i.e., 
77 = 0. 

The magnetic torque makes the disk unstable, while mag- 
netic diffusion has a stabilizing effect, suppressing instability 
at low f2 and/or torque (see Figs. 3 and 4). Comparing Figs. 6 
and 9, we see that the disk is less unstable in low temperature 
case (H = 0.001) while the instability occurs for lower values 

of Kq. 

To see the effect of magnetic diffusion, Figs. 10 and 11 
show the growth rates for the cases with fj = 0. In the absence 
of magnetic diffusion, the disk is always unstable, though the 
growth rates can be very small. 

The physical reason for the instability found here has been 
described in LPP (see discussion in section 2). A perturba- 
tion increasing v r causes the poloidal field to be bent close 
to the disk plane, as the field is advected inwards by the ac- 
creting matter. This tends to increase the mass flow along the 
field. There is also an opposing effect, however. As the poloidal 
magnetic field is bent towards the disk plane, the radial cur- 
vature force on the disk (opposite to gravitational force in- 
creases). Due to this increased support against gravity, the ro- 
tation rate of the field line decreases. As a result, the "po- 
tential barrier" for mass flowing along the field line is larger, 
decreasing the mass flux and magnetic torque. The instability 
decreases with increasing field strength through its effect on 
the rotation rate, since the instability is driven by the magnetic 
torques. We find that the inclination of the magnetic field still 
has the strongest influence on the magnetic torque, and domi- 
nates over the change in the potential barrier. 

Since the instability is caused by the magnetic torque by 
the centrifugally driven wind from the disk, the instability dis- 
appears as the magnetic torque vanishes (Figs. 3 and 4). The 
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Fig. 11. Same as Fig. 8, but without magnetic diffusion, fj = 0. 

disk is stable without a magnetic torque, all three modes are 
then damped by magnetic diffusion. If the magnetic torque is 
large, the instability time scale is comparable with the dynam- 
ical time scale of the disk, and the instability is not much af- 
fected by magnetic diffusion (compare Figs. 6 and 10). We find 
that the instability is insensitive to the perturbation wavelength. 

Associated with the perturbations in the disk are changes of 
the field configuration above the disk. These affect the acceler- 
ation of the wind, and hence the wind torque. In our analysis, 
we have taken the response of the wind torque to these changes 
to be instantaneous. This is motivated by the fact that the mag- 
netic torque acts on the disk surface, and changes as soon as the 
azimuthal field component at the surface changes. This is per- 
haps somewhat contrary to the impression that one might have 
based on the steady wind model. In the steady wind model, the 
torque is often pictured as 'effectively acting' at the Alfven sur- 
face. In reality, the torque is constant along the field line, hence 
acts equally at the disk surface. When changes due to motions 
in the disk take place, changes in the torque travel up at the 
Alfven speed and are first felt in disk itself. 

The actual propagation of these torque changes has not 
been considered here, since its time scale is short, of the or- 
der of the Alfven travel time over the disk thickness. The time 
scale TAd for an Alfven wave to travel over the scale-height 
is of the order TAd ~ c s /Va^ _1 , which is short except for 
weak fields. The time scale for Alfven wave traveling from the 
disk surface to the Alfven point in the wind is in fact shorter 
than this, ~ l/fi - The disk- wind coupling time scale is there- 
fore of the order of the disk dynamical time (Wardle & Konigl 
1993; Konigl & Wardle 1996). The present analyses are valid if 
the instability growth time scale is much longer than the disk- 



wind coupling time scale. This is the case since for weak fields, 
where the coupling time is long, the growth time of the insta- 
bility itself is also long. 

In order to check on the effect of finite Alfven travel times 
in the disk-wind coupling process on the instability properties 
of the disk, we have done a test calculation with a simplified 
description of this effect. In this calculation, we manually in- 
duce a time delay of the dynamical time scale of the disk on the 
wind torque: 

dT m (t) = T mQ -T m (t) ^ m 

9t Tdclay 

where T m0 is given by Eq. (15), Td c i a y is the assumed disk- wind 
coupling time scale. We then add this equation to the set of disk 
equations. We find the growth rate of the instability becomes 
slightly lower than before, but the instability properties have 
not been altered qualitatively. 

The nonlinear evolution of the instability analyzed gas to 
be studied by numerical simulations. Some such simulations 
have already been reported by Agapitou & Papaloizou (2001, 
in preparation). 
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